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SUMMARY 

The theory and conditions for stability of a satellite 
with a 24- hour period are given. Bohlin’s resonance theory 
was applied to obtain the solution. It is shown that the in- 
tegrals of the problem can be represented in series form, 
with respect to the small parameter w, which is propor- 
tional to the mean motion of the critical argument in a non- 
resonance case. Expressions for the period of libration and 
the mean motion of the critical argument in the unstable 
case are also given. A system of formulas is presented 
which can be used to compute any particular case. 
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ON THE MOTION OF A 24-HOUR SATELLITE* 


by 

Peter Musen and Ann E. Bailie 
Goddard Space Flight Center 


INTRODUCTION 

In this paper the authors have investigated the stability conditions of a satellite with a period of 
revolution approximately equal to one day. The criteria for stability are obtained in a form valid for 
large inclinations- provided that the eighth power of the eccentricity is negligible. 


THE DISTURBING FUNCTION 

The disturbing function consists of a secular part produced by the zonal harmonics k 2 and k 4 ,and 
a periodic part produced by the ellipticity of the earth's equator. The periodic disturbing function has 
the form 


3 ^ A 22 / - 2 x 
F* = ~“-7- (x 2 " y 2 ) 


1 2 r 5 


(1) 


The x-axis is directed along the semi- major axis of the equator, and the z-axis is directed along the 
axis of rotation of the earth. Substituting 

x _ 1 + 2 ° S ~ cos (f+w + Q- n't) + 2° S 1 cos + * 


y 1 + cos i . , „ _ . . 

— = 2 sin ( f + Co + Q * n t ) “ 


1 - cos i 


sin(f + to - 0 + n ' t ) 


sin{f+a)) sini 


•This report has been published in substantially the same form in ]. Geophys. Res. 67(3):1123-32, March 1962. 
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into Equation 1, we obtain 




3 ^ A 2 2 

8 ^ 


( 1 + cos i ) 



cos ( 2f + 2 <*j + 2ft ~ 2n ' t } 


3 ^ A 2 2 

+ 


i cos (20- 2n' t ] 


3 M A 22 / a \ ^ 

'g — (1-cosi) 2 (“j cos ( 2f + 2oo - 20 + 2n ' 1 ) 


( 2 ) 


Developing Equation 2 into a series in terms of the mean anomaly l with coefficients developed in 
powers of the eccentricity, and retaining the long period terms only, we deduce 

F j = Q 0 cos (2 1 + 2co + 20 - 2n ' t } 

+ Qj cos ( 2l + 20 - 2n r t ) + Q 2 cos ( 2 1 ~ 2co + 20 - 2n ' t ) * /g\ 


where 


3 ^22 , ,, ( 5 . 13 . 35 ,\ 

Q 0 - 5 ~^r ( 1 + cos i ) 2 ^1 - 2 e 2 + 16 e ~ 288 e J 


3 ^22 

Qi - 4 ~r s 


,.(9 , 7 , 141 

m 1(4 e + 4 e + 64 e j 


3 ^22 , ./ 1 , 7 6 \ 

Q 2 - g--^~ (1-cosx) [24 e + 240 e j 


If the mean motion of the satellite is such that it causes the satellite to remain above a particular 
longitude of the earth for some time, the first term in Equation 3 is the most significant, and will be 
treated in accordance with the theory of resonance. The last two terms will produce only small, long 
period terms in the osculating elements. 

The coefficient A 22 is approximately of the same order as the coefficient of the fourth zonal har- 
monic, and consequently, the secular part F 0 of the disturbing function need not be developed beyond 
the results established by Brouwer (Reference 1). Adding the term n ' y^ta( 1-e 2 ) cos i, produced by the 
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rotation of the earth to Brouwer's development, we have 


M j ^ 4k 2 / 1 3 H 

F n = + n 11 + 

0 2L 3 


L 4 G 


3,3 I 2 + 2 r 2 


^4/15 1/ 2 — - to — + — — 

T 15 ” V 16 G 7 16 G* A G 2 3 G 4 


M 6 k 2 2 


10 


15 L 5 / 18 H 2 

+ 33 1 - T G 2 + G 4 


3 L° 


H 2 _ H 4 


+ 8^ U-6 ^ 4 ^ 


G 4 


15 L 7 ( H 2 H 4 
32 l 1 " 2 G 2 ' 7 


where L, G, and H are Delaunay variables 


1 cos 1 , 


L = y^Ii , G = Aa(l- e 2 ) , H = ^ /xa 1 ~e 2 ) ■ 

and k 2 and k 4 are the coefficients of the second and fourth zonal harmonics, respectively. 

Considering the form of the main argument, it is more convenient to use the canonical set of 
Poincare rather than that of Delaunay. 


Thus 


yjja 


yTZS (l - y 1 - e 2 ) = L - G, 


and 


= |/7Ia (l ~ }/l “ e 2 cos i ) = L - H, 


- l + co + 0 “ n ' t , 


y 2 - - a;, and 


y = - n + n't. 


A 


> (4) 


The complete Hamiltonian is 


F = F 0 + Fj , 


(5) 


but the disturbing function can be expressed more conveniently by means of the auxiliary quantities 
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x r , e, and 7, where e - x 2 /x l , y = l-x 3 y'x 1 . And the following relations exist: 


L * Xj , 


L = 1 * * * 


H _ . _ 7 

~Q - COS l ~ Y - 


e 2 = 2 e - € 2 . 


The following expressions were used in the actual computations: 


? 0 ' n y*\ 


M k 2 

+ (l + 3e + 6 e + 10 e 3 ) + y 2 (3 + 15e + 45e 2 + 105e 3 )J 


M Kj 3 3 

+ Yo~ 52 (4 + 14e + 19e 2 - 21e 3 ) - jg y 2 (l 6 + 114e + 459e 2 + 1371e 3 ) 


32 y 4 (76 + 790e + 4515e 2 + 18755e 3 ) 


M 6 k r 

^ (2 + 20e + 95e 2 + 315e 3 ) - y 2 (2 + 24e + 141e 2 + 573e 3 ) 

! L 


+ yg y 4 (l4 + 196e + 1365e 2 + 6545« 3 ) 


^ a 22 


gg (36 - 180e + 207e 2 - 152e 3 ) +^ 7 (36 - 144e + 63e 2 - 89e 3 ) 


^ y 2 (36 - 108e ~45e 2 - 134e 3 ) 


m 4 A 22 r 3 3 1 

r~ 32 i 366 + 38e2 + S5€ 3 )-- 32 7 2 (36e + 110e 2 + 269e 3 ) * 

Xj L J 

/^ 4 A j - 

|SS ( S£j + 2 « 3 ) - ^ y (s « 2 + 7e 3 ) + jjjjj y J (5e 2 + 12e 3 ) 
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The first and second derivatives of F 0 and Q 0 with respect to x, can be obtained by using the expression 

<9F 1 - y dF £ d¥ 

F ' = 3T, + V Sy - ^ JI ' 

from which we obtain: 

F ' = n > - + — .1 [(2 + 7e + 16e 2 + 30e 3 ) + 2y (l + 5e + 15£ 2 + 35e 3 ) - y 2 {s + 45e + 150e 2 + 385£ 3 )] 

0 x. 3 2x. 7 


^ TT [33 (40 + 154e + 228e 2 - 273e 3 ) + |- y ( 16 + 114e + 459e 2 + 1371e 3 ) 

- ^ y 2 (64 + 494e + 2142e 2 + 6855e 3 ) “ J y 3 (76 + 790e + 4S15e 2 + 18755e 3 ) 

+ y 4 (1064 + 11850e + 72240e 2 + 318835e 3 ) 


[|| ( 4 + 44e + 228e 2 + 819e 3 ) + “ y (2 + 24£ + 141e 2 + 573£ 3 ) 


- ^ y 2 ( 8 + 104e + 658f 2 + 2865e 3 ) - ^ y 3 (2 + 28e + 195e 2 + 935e 3 ) 


9 P 

+ Ye r 4 (28 + 420e + 3120e 2 + 15895£ 3 ) 


&L [3 ( 6 + 23e + 57c 2 + 115e 3 ) + 6 y (8 + 45e + 1S0 £ 2 + 38Se 3 ) - 9y 2 (l2 + 75e + 275e 2 + 770e 3 ) 

x, 4 x, 8 L 


F o" 


(l88 + 696e + 564e 2 - 4653« 3 ) + |y (l 28 + 988 e + 4284e 2 + 13710e 3 ) 

1 

- | y 2 ( 340 + 2668e + 11550£ 2 + 36085e 3 ) - -f y 3 ( 1064 + 11850e + 72240e 3 + 318835« 3 ) 

+ || y 4 ( 532 + 6320f + 409366 2 + 191301e 3 )J 


^ IT -T ( 6 + 72f + 400e 2 + 1529£ 3 ) + 4| y (8 + 1046 + 658f 2 + 2865£ 3 ) 


- ^ y 2 ( 142 + 1988f + 13440e 2 +622156 3 ) - ^ y 3 (28 + 4206 + 3120 e 2 + 158956 3 ) 


+ ^ y 4 ( 14 + 224£ + 1768£ 2 + 9537e 3 ) 
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M 4 A, 


43 (72 - 486e + 765e 2 - 595e 3 ) + -| 7 ( 6 - 29e + 17e 2 - 21« 3 ) 




m 4 a. 


48 


+ 53 y 2 (l44 - 486e - 225e 2 - 737« 3 ) 

(288 - 2844e + 6273« 2 - 5194e 3 ) + ^y (720 - 4212e + 3285e 2 - 3421s 3 ) 

+ 13 7 2 ( 432 - 1620e - 825s 2 - 2948e 3 ) 


If only the secular and the first, most important, periodic terms of the disturbing function are 
retained, we have: 


F - F o + % cos 2yj . ( 6 

The libration points and the points lying on the intersection of two branches of the separatrix are 
determined from the equations 


8F _ aF o d % 

= dx l + cos 2 y\ 


and 


dy x ~ ~ 2 % sin 


if the values of x 2 and x 3 are fixed. From Equation 6 we deduce that 


3F 0 

3Q 0 


3*1 

+ Bxj ~ ® 

(7) 


for y, - 

0, 77 , 

(8) 

Bxj ” Bxj 

= 0 

(9) 

for y x - 

77 377 


~2 ’ ~2 

(10) 



The values y = 77/2 , 377/2 reduce the disturbing function (Equation 5) to a minimum. Consequently, 
they correspond to the stability position and determine the libration points. The set y x = 0,77 gives 
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the points on the separatrix and corresponds to unstable positions. For satellites moving in the 
equatorial plane the stable positions are on the equator's minor axis, and the unstable positions are 
on the major axis (Reference 2). It is convenient in the complete problem, as defined by Equation 5, 
to retain the development around the characteristic points defined by Equations 7 through 10, and to 
write the disturbing function in the form 


F = R„ + R, 


where 


R 0 = F o - Qo ’ 

R, = 2Q 0 cos 2 yi + Q, cos (2y t + 2y 2 ) + Q 2 cos (2y, +4y 2 ) 

for the stable case, and 

R 0 = F o + Qo ’ (11) 

R, = - 2Q 0 sin 2 yj + Q, cos (2y t + 2y 2 ) + Q 2 cos (2 yi +4y 2 ) ^ 

for the unstable case. 

The canonical transformation (Equation 4) removes the time and the argument y 3 from the dis- 
turbing function. Consequently, this problem contains the energy integral 


R„ + R, 


- C 


(13) 


and the integral 


= yrn ^1 - - e 2 cos i ^ 


is constant. 

THE STABLE CASE 

By substituting x 2 = as/a yi , and x 2 = ds/dy 2 into Equation 13, we transform it into a Hamilton- 
Jacobi partial differential equation. Letting 

s = s 0 + s, (y x ) + s 2 ( yi , y 2 ) + s 3 (y, , y 2 ) > 
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+ Cj + C 2 + C 3 + 


and 


s o = 


+ a 2 y 2 + a 3 y 3 


where a, and a 2 can be considered as two constants' of integration. We replace the partial differential 
Equation 13 by the system of Bohlin's equations (Reference 3). These may be deduced by developing 
the left-hand side of Equation 13 into a Taylor's series about a 3 , a 2 , a 3 : 


R 0 (a,, <iji a 3 ) - - C 0 , 


(14) 


, , <5S. V 

a + 2 w Jo\57;/ + w 10 ^ + 2Q 0 cos 2 y, + Q,cos (2yj + 2y 2 ) + Q 2 cos (2y, + 4y 2 ) = 0- (15) 


as 2 ! 

w oi dy 2 + 2 W 20 \dy x 


dS % 


\ 9S 2 <9 S j 


i f* s i 


4 ... 


01 + V w ‘° +W20 ^7 ) W\ + 3y, dy 2 w »» + 6 w »\3y„ 

[ 2Q 0 ' cos 3 yj + Q,' cos (2y t + 2y 2 ) + Q 2 ' cos (2y t + 4y 2 )J = - C, > (16) 


3 S f 


where 


a ltjR 0 (a,, q 2 , a 3 ) 

b da j 1 daj 

and 

JQi ( a l< “ 3 ) 

^1 da t 

Equation 14 can be considered as a defining equation for C 0 . 

Imposing on our solution an additional condition — that no secular term with respect to y 2 is con- 
tained in S„ S 2 , S 3 , , but only in S 0 -we deduce from Equations 14, 15, and 16 that 


Si = *, (y,) - 


^2 (Vl) 


_Qi_ 

2w nl 


sin (2y, + 2y 2 ) - 4W(m sin (2y, + 4y 2 ) 
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and also 


T w 20 $1 2 + w » $1 + 2Q 0 cos 2 y, - 0 • (17) 

*2 ( W 10 +W 20 *!*) + i W 30 *1 3 + Qo' COs2 yi = " C 1 • ( 18 ) 

From Equation 17, we have the standard first approximation 

4>i = “ w + A . (19) 


where 


W 


W 

w 


10 

20 


and 



4Qo 


w cos y , 
W 20 


(20) 


If 


4Qo 


> 1 


( 21 ) 


then cosy, oscillates between the limits "( w / 2 ) (y' w 2 o/Oo)- + ( w / 2 ) (y/ w 2 </Oo) • This condition (Equation 21) 
can be written in the form: 


(F 0 '-<V) 2 - 4Q 0 (f 0 ”-q 0 “)< 0 - (22) 

(where, in this case, the primes represent differentiation with respect to a, ), together with the 
condition 


. (V-Qq ') 2 

cos y » < 4Q 0 (F 0 " - Q 0 ") 


which, for a certain moment of time, must be fulfilled for the motion to be stable. By eliminating <j> l 
and cos J y, from Equation 18 by means of Equations 19 and 20, we can write an expression for <f> 2 ‘ : 




( 1 w 30 1_ 

\6 w 2 o 2 Q 0 


A 2 + 


'30 


Qo' 

Q 0 


Aw - 


1 

2 




1 W 30 
6 W 20 


+ 
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The divisor A in the last term may 
continuity, we set 


10 


( 30 ) 



we obtain the Hamiltonian function S with its integrals in the form: 


and 


S = (aj+AOJy, +a 2 y 2 + a 3 y 3 ± J ^ A«> - A< 3 > cos 2 y , dy , 

+ A (4) sin 2 y j + A (5) sin ( 2 y j + 2 y 2 ) + sin ( 2 y x + 4 y 2 ) 


x - = (a,+A^)) ± |/a ( 2) - A< 3 > cos 2 y t 




as 

‘a " ay 2 


+ 2 A< 4 > cos 2 y x + 2 A (5) cos ( 2 y, + 2 y 2 ) + 2 A ( ^ cos ( 2 y x + 4 y 2 ) 


= a 2 + 2 A CS > cos ( 2 y t + 2 y 2 ) + 4 A <9 cos ( 2 y, + 4 y 2 ) 


as 

x 3 “ a y 3 " a 3 


as 

da . 


n j t + / 3 j = (l+A/») y j ± Mj 


! 


dy, 


ypjfi - cos 2 yj 


± \ - A< 3 > cos 2 y, dy 3 

+ A,^sin 2 y, + A, (5 > sin ( 2 y , + 2 y 2 ) + A, (6 > sin ( 2 y , + 4 y 2 ) 1 


§§- n 2 1 + = A J (1> yi + y 2 1 M 


j 


dy, 


/a<* - A (3) cos 2 y j 

1 a w r / 

± "J J / A (2) - A<^> cos 2 y , dy , 

+ A 2 < 4 >sin 2 y, + A 2 (5 > sin ( 2 y, + 2 y 2 ) + A 2 (6) sin ( 2 y, + 4 y 2 ) ’ 


dS 
da - 


n 3 4 + = A 3 (,)y i + y 3 1 m 3 

A < 3 > 


i- 


dy, 


yv 2 * - a < 3) cos 2 y j 

1 T J VA<^ - A* 3 * cos 2 y , dy , 

A 3 W sin 2 y, + A® sin ( 2 y, + 2 y 2 ) + A 3 <6 > sin ( 2 y, + 4 y 2 ) 


(31) 


(32) 


(33) 


(34) 


(35) 
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In this representation the terms containing y 3 in the argument are of the first order with respect 
to the paramater k 2 . In addition, these terms contain even powers of the eccentricity as a factor and 
consequently, will be small from the start; the computation of the second order term in k 2 with the 
argument y 2 shows that it can be neglected. All quantities can be considered as functions of n, t +/3., 
(i = 1,2,3) but, of course, only y 2 and y 3 will have a secular term, since yj does not possess any 
such term. 

It is not difficult to continue the process of computing s, if necessary. However, taking present 
day knowledge of the numerical values of geodetic parameters into consideration, it was found that 
even the development of Equations 31 through 35 proved to be accurate - from the practical point of 
view -overly accurate. In this solution R a was originally considered to be of the second order with 
respect to w 10 . It must be pointed out that this classification is purely formalistic, and loses its 
significance after the development is completed. The important characteristic of the solution is that 
it can be developed into a series in w, with the coefficients depending upon a 2 , a 2 , a 3 . The develop- 
ment is not made in powers of^lTJ as might be expected. This feature was observed initially bylzsak 
(Reference 3) in his solution of the critical inclination problem. 

The method presented here does not introduce the small divisor A in the determination of <P n which 
appears in the expression for S n : 

s n = <p n (y,) + trigonometric terms in y, and y 2 . 

Every <£ n ' will be a polynomial in A, if the constant of energy C is decomposed properly, to remove 
the poles with respect to A. This can be easily shown by applying the "from n to n + l proof," since 
the equation for the determination of <p n ' has the form: 


A + P 




<P 1 1 ; 


ly i) 


= c - 1 


where P n is a polynomial in ,<£ 2 ' , * • • , <P n f . x cos 2 y 1( The elimination ofcos 2 y 15 by means of 
Equation 17 and the proper determination of C n ^ ,will lead to the representation of <p n ' in polynomial 
form, providing it has been shown that <P X , <£ 2 ' , * • * , <P n L x are polynomials inA # Eliminating higher 
powers of A in favor of cos 2 y 1 , we deduce that 


<*>„' = a n (c°s 3 yj + A/3 n (cos 2 y j) • 

where a n and are polynomials in cos 2 y 2 with polynomial coefficients inw. This result is similar to 
that obtained by Izsak (Reference 4) for the critical inclination problem. For the partial derivative 
of S n with respect to y 2 we have 

dS 

ij; = T n (D) (y»- y») + AT ^ ( y i- y a) - 

where T n <0) and T n (1 > are trigonometric polynomials in 2y, and 2y 2 with polynomial coefficients inw. 
When integrated, this has been found to be purely trigonometrical with respect to y 2 . 
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In the stable case y x will have a long period oscillation about 77/2 or 3 tt/ 2. It is also of interest 
to know the approximate period of this deviation in longitude. From Equations 23 through 29, we have 


A/ l > 


1 - i £ • ♦ ° M 


(36) 


A,< 2 > 


= 2w + 




0 (w 3 ) > 


(37) 


A/ 3 > 

"a® 


0 (w) ' 


(38) 


M, 


1 

w + 2 



+ O (w 3 ) ■ 


(39) 





(40) 


Neglecting small long-period terms, we can write, from Equation 32, the expression for the period of 
libration T : 


where 


Setting 


n,T = 2(l + A 1 (,) ) arc cos (•£■) “ arc cos (-■£) 


M, dy t 

^j«.-VlA) + k2 C ° S2 y i 


A< 3 > 




AW 


/■cos ( 1A ) 

V- 

J cos (-1/k) 


1 - k 2 cos 2 y j dy j » 


a (3 > _ 4 Qq 

A (Z) w 2 w, n 


k cos y j - sin u , 
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we have 


-cos 1/k ) 


^cos _1 (-l/k) 


d yi 

j/l “ k 2 cos 2 y j 


/•TT/2 



and 


y.cos“ I (l/k) 

•7T/2 

j/l - k 2 cos 2 y t dy , = - £ 

cos 2 udu 2 „ fl\ . T fl\ fl\ 

, x = ’ k K IkJ " 2k [e U-k (JJ 

Jcoi _1 (-l/k) J 

o y 1_ ^ sin2u 


where K(1A) and E(iA) are the standard elliptic integrals of the first and second kinds having 
the modulus 1/k. In terms of the new variables, 


n i T 


= 2 


1 + 2 arc cos (j~j - tt ~ y 


1 A « 


Yam 




(I) 


A,® 

2k 


A<3) 


E (r) ~ K (y) 


(41) 


Substituting Equations 36 through 40 into Equation 41 yields: 


2 W 30 

n , T = -3^w 


2 cos” 


( 4 ) - -] - 4 (i K ® + o (t) 


Equations 14 and 23 give the expression for c = C 0 +C X , from which 


ac 

1 “ 3a x 


i . _J_ (?» 
2W 20 \ W 20 



w 


and 



In order to find the period in the vicinity of the libration point, we expand Equation 42 in powers 
of 1/k. Since w approaches zero as the libration becomes smaller, 
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and for small values of 1A 


.Jl) . s.i.U 

\k/ 2 k 6 1,3 


/ 1 \ 7 7 / 11 

Ik) “ 2 ( 1 + 4 k 2 


l +x nr + 


Consequently we deduce that 


2 w 30 Q o' 1 /So Qo” 

^ w 3 TT W 20 Q 0 2w 20 l w J0 Qo 


UNSTABLE MOTION 

If the condition (Equation 22) is not satisfied, the motion is unstable and the Hamiltonian function 
takes the form 

F = R 0 + R, - 

with R 0 and defined by Equations 11 and 12. Performing operations similar to the previous ones, 
we find that 

S = (a, + B<‘)) y, + a 2 y 2 + ^ y 3 + J j/b<* + B® sin 2 y, dy 2 

+ B< 4 > sin 2y, + B <5) sin (2y 4 + 2y 2 ) + B (6) sin (2y 2 + 4y 2 ) • 

where 


... 2 , w 30 

B (l) = -W-TW 1 

J 20 


1 W 30 %' \ Q 0 


W 20 Qo / W 20 


- w 2 1 + 


1 l W 30 Qo' 


2 \ w 20 Qo 


B<3) = 


4Qo f 1 / w 30 Qo' , \l 2 


W 20 1 + 2 \ W 20 W, 


w >0 


_ / 1 ^30 _ 1 Qp'\ Qo 
\ 7 W 20 2 Qq / W 2 0 
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and 


B<s> = 


Qi 


2w 


01 


B(«) = - 


4w 


01 


Despite the similarity between some A f s and B T s, they are not identical since different values of 
R 0 are used in each case. The constant of energy has the same analytical form as in the stable case: 


( 1 W 30 1 ^0 \ 

c " R o (V a i) + U^ + 2 id* 3 


( 46 ) 


Thus we have, as in the previous case, 


- a 1 + + ^/bC 2 ) + sin 2 y ~ + 2B W cos 2y x + cos (2yj + 2y 2 ) + 2B^ cos (2y t + 4y 2 ) * 


x 2 - a 2 + 2B^ cos (2y, +2y 2 ) + 4B<$> cos ( 2y j + 4y 2 ) 


and, putting 


B/ 1 ) 


<iB« 


dC 
da . 


, , B.» 

N = •*• B .«> -B« - J — 

I 2 \ 1 g(3) 


we have, designating by the additive constants of integration /?,, , p 3 


n t t + /3, = (1 +B 1 ^)y 1 + 


"'1 


dy, 




1 B 1 (3) f n — 

f=^ + 2 B (3)J' + B(?>sin2 yi dy, 


sxn-y, 


+ B,W s in 2y, + B, <s) sin (2y, + 2y 2 ) 


+ B,W sin ( 2y, + 4y 2 ) ; 


( 47 ) 
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n 2 t + P 2 = b 2 « y t + y 2 + 


4' 


d y t 

|/b^ 2 ) + B® sin 2 yj 


1 ^ f f 

+ y 


Q( 2 ) + sin 2 yj dyj 


+ B 2 < 4 > sin 2y j + sin (2y x + 2y 2 ) + B 2 (6 > sin (2y t + 4y 2 ) * 


n 3 t + p z ~ B 3 (1 > y 1 + y 3 + N 


f dy i B W r 

I — ~ + -x — ~ i/b® + B< 3 > sin 2 y dy 

j B® + B^ 3 > sin 2 y j 2 


+ B« sin 2y, + B 3 < 5 > sin (2yj + 2y 2 ) + B® sin (2y t + 4y 2 ) 


In the unstable case 


B( 2 ) + B^ sin 2 y 1 > 0 


(48) 


and the square root of Equation 47 (and its reciprocal) may be developed into a Fourier series in 
2y 3 ; the argument yjWill possess a real secular term, which is absent in the stable case. The 
coefficient of y, in the right side of Equation 46 is 


P = 1 + B + 


2Nj n/2 

7 T 

0 

J 


dYi 


3,d) ,tt/2 


+ B <3 > sin 2 y j 


1 r n/ 1 t 

J ya» + sin 2 Yj dy, 


(49) 


Substituting in the integrand 


yi 



and 


. a _ 

B (2 > + B& 

we reduce the integrals to the normal form and Equation 49 becomes: 


P = 


2KN. E 

+ B,d) + 

„■/&*> + B< 3 > 77 


B.d) , 

V Bd) + g(?) , 

B 0) ' 


(50) 
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where 


K(k) = 


r 7T/2 d<£ 

J o i/l - k 2 sin 3 


E( k } = x/l — k 2 sin 2 0 d^> 


Taking Equations 43, 44, and 45 into account, we further deduce that 


B ro = _ ! _ i w _J1 + 0( 'w 2 ) , 




= 0(w) 


1 / w 30 ^0 \ / 3 \ 

N > = w + 2 \ ^ -3 Qr) w + °( w ) 


It follows from Equations 50 through 54, neglecting the terms of the second order in w, that: 


1 W 30 2k' K / Qo' 


p = - 3 + — u- TtJ* 


rfhere k' is the complementary modulus, k 2 + k' J = l, and 


l/ W 30 <?o' 


t/ B ( 2> + g(3) 1 + 2 \ w JO Qoj* 


From Equation 46, neglecting terms of the third order, we have 
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Using Equations 54 and 55, we can compute the mean motion v, = n,/P of the argument y, with an 
accuracy up to the terms of the second order in w. The value of denotes the speed with which the 
satellite will depart from its original position over the earth's surface during the course of time. 

CONCLUSION 

The theory of motion of a 24-hour satellite under the influence of the ellipticity of earth's 
equator has been developed using a resonance theory. The expressions for the elements of motion 
can be represented in the form of a series with respect to the parameter w , which would be closely 
associated with the mean motion of the main critical argument (n -n') t +n + <a in a nonresonance case. 
Canonical elements and the Hamilton-Jacobi partial differential equations were used to solve the 
problem. This method of solution was chosen because of its flexibility with respect to the form of the 
integration constants which are adjusted so as to remove small divisors from the solution. 

The conditions are established for stable and unstable types of motion; no severe restrictions are 
imposed on the values of inclination or eccentricity. The formulas are developed to the point where a 
numerical development can be easily accomplished for any particular case. Using this method, an 
extension of Hori’s critical inclination theory (Reference 5) can be easily obtained. The method 
described herein can also be applied to a more general case. If the ratio of mean motions is 
n/n' = 1/p, and not n/n' = 1 as in our case, the theory given here can be easily extended by using the 
set of canonical variables: 

x = , y ^ = p l 0 + O) — n't * 

X 2 = - e2 ) ’ y 2 = ~. a> ' 
x 3 = - y'l-e^cos i) • y 3 = - n + n't - 

Comparing this theory with observations can help to better determine the coefficient of ellipticity 
of the earth's equator. 
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